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ABSTRACT 

Context. Interstellar filaments are an important part of the star formation process. In order to understand the structure and formation 
of filaments, the filament cross-section profiles are often fitted with the so-called Plummer profile function. Currently this profiling is 
often approached with submillimetre studies, especially with Herschel. If these data are not available, it would be more convenient if 
filament properties could be studied using groundbased near-infrared (NIR) observations. 

Aims. We compare the filament profiles obtained by NIR extinction and submillimetre observations to find out if reliable profiles can 
be derived using NIR observations. 

Methods. We use J-, H-, and K-band data of a filament north of TMC-1 to derive an extinction map from colour excesses of 
background stars. We also use 2MASS data of this and another filament in TMC-1. We compare the Plummer profiles obtained 
from these extinction maps with Herschel dust emission maps. We present two new methods to estimate profiles from NIR data: 
Plummer profile fits to median A v of stars within certain offset or directly to the A v of individual stars. We compare these methods by 
simulations. 

Results. In simulations the extinction maps and the new methods give correct results to within ~ 10-20% for modest densities (p c = 
10 4 -10 5 cm -3 ). The direct fit to data on individual stars usually gives more accurate results than the extinction map, and can work in 
higher density. In the profile fits to real observations, the values of Plummer parameters are generally similar to within a factor of ~2 
(up to a factor of ~5). Although the parameter values can vary significantly, the estimates of filament mass usually remain accurate to 
within some tens of per cent. Our results for TMC-1 are in good agreement with earlier results obtained with SCUBA and ISO. High 
resolution NIR data give more details, but 2MASS data can be used to estimate approximate profiles. 

Conclusions. NIR extinction maps can be used as an alternative to submm observations to profile filaments. Direct fits of stars can 
also be a valuable tool in profiling. However, the Plummer profile parameters are not always well constrained, and caution should be 
taken when making the fits and interpreting the results. In the evaluation of the Plummer parameters, one can also make use of the 
independence of the dust emission and NIR data and the difference in the shapes of the associated confidence regions. 

Key words. ISM: Structure - ISM: Clouds - Stars: formation 



1. Introduction 

Interstellar filaments have long been recognised as an impor- 
tant part of the star formation process. Recent studies have 
confirmed that filamentary structures are a common feature in 
dense interstellar clouds and different stages of potentially form- 
ing stars are located in these filaments (e.g. | Andre etaLI 1201 Ot 
lArzoumanian et all 1201 it [Juvela et aLll2012bl) . Therefore, it is 
important to study the formation and structure of the filaments 
in order to understand the details of star formation. 

The mass structure of molecular clouds can be studied via a 
number of methods. These include molecular line mapping, ob- 
servations of dust emission at far-infrared/submillimetre wave- 
lengths, star-counts in the optical and near-infrared wavelengths, 
and measurements of colour excesses of background stars to 
produce an extinction map. The surface brightness based on 
scattered NIR light has also proved to be a potential method 
for studying mass distributions (see e.g. iPadoan et al.l [2006; 

* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



lJuvela etldl 120061 120081) . However, all techniq ues have their 
own drawbacks (see, e.g.. iGoodman et aT1l2009l for a compari- 
son of several methods). For example, line and continuum emis- 
sion maps are subject to abundance variations (gas and dust, re- 
spectively) and variations in the physical conditions. Mass esti- 
mates based on dust emission and on the estimation of colour 
temperature can also be biased because of line-of-sight temper- 
ature variati ons, especially in potentially star forming high den- 
sity clouds dMalinen et alJl201 ll) . 

The details of filaments have been studied for a long time 
using various metho ds, including extinction maps (see e.g. 
ISchmalzl et al„ 2010, and references therein). However, the cur- 
rent interest in studying fi lament profiles originates from submil- 
limetre observations (e .g.lNutter et al.ll2.Q08l; lArzoumanian et al.l 
1201 lb iHill et al.ll20Tll lJuvela et al.ll2012bl) . espec ially the high 
resolu tion data of Herschel Space Observatory dPilbratt et al.l 
120101) . Many recent studies have fitted the filament cross-section 
profiles with the so-called Plummer profile function (see details 
in Section|4]i, because the parameters of the function can give in- 
formation about the formation and structure of the filament (see 
more details e.g. inlArzoumani an et al ,l201ll:lJuvela et al.l2012al 
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submitted). Large magneto-hydrodynamic (MHD) simulations 
have shown that filaments ca n be formed as a n atural result 
of inter stellar turbulence (e.g. iPadoan et al.ll2001l) . lJuvela et al] 
d2012al submitted) have studied the observational effects in de- 
termining filament profiles from synthetic submillimetre obser- 
vations. 

The use of far-infrared/submillimetre wavelengths is more 
expensive, as spaceborne satellites are needed. However, in the 
future groundbased large-area surveys made with submillime- 
tre wavelengths ( for example with SCUBA-II) will also become 
available (see e.g. lWard-Thom pson et al. 2007). If these data are 
not available, it would be more convenient if groundbased NIR 
observations could be used to study the detailed properties of 
cloud filaments. The use of dust extinction could also be more 
reliable, as it is not dependent on the easily biased estimates 
of colour temperature. It would also be useful to compare mass 
structures obtained using different methods. We therefore aim to 
compare the filament Plummer profiles obtained from NIR ex- 
tinction data to those obtained from submillimetre dust emission 
observations. 

The Taurus molecular cloud is one of the closest (dis- 
tance ~ 140 pc) and most studied star-forming regions (see e.g. 
iNutter et al.1 12001 ISchmalzl et al.ll20Tot iLombardi et al]|2010l 
and references therein). The cloud has also been mapped with 
He rschel as part of t he Gould Belt Survey (iAndre et al.ll2010l) . 
see lKirk et alj d2012f) and Palmeirim et al. (in prep.). It is there- 
fore a good candidate for studying the structure of filaments. 
In this paper, we will examine in detail two cloud filaments in 
Taurus, comparing the results obtained from NIR and submil- 
limetre data. Additionally, we examine two more distant clouds 
mapped by Herschel. 

The contents of this article are as follows: Observations and 
data reduction are described in Sect. |2]and the methods to derive 
column densities in Sect. [3] We present our methods to derive 
filament profiles in Sect. [4] The results of profiling observed fil- 
aments are shown in Sect. [5] We examine the performance of 
the analysis methods in Sect. [6j with the help of simulations. 
We discuss the methods and results in Sect. [7] and present our 
conclusions in Sect. [8] 

2. Observations and data reduction 

We study two filaments in t he Taurus molecul ar cloud, TMC- 
1 (also known as Bull's tail (Nut ter et al.l l2008)) and a filament 
north of TMC-1, which we call TMC-1N. The central coordi- 
nates are RA (J2000) 4h41m30s and Dec (J2000) +25°45'0" for 
TMC- 1 and RA (J2000) 4h39m36s and Dec (J2000) +26°39'32" 
for TMC- IN. A combined extinc tion map o f the tw o fields is 
shown in Fig. Q] According to iRebull et al.1 (1201 ll) . there are 
no known young stellar objects (YSOs) in TMC-1N. There are 
several YSOs in TMC-1, however, six of them in the filament 
analysed here. YSOs are redder than background stars and could 
therefore bias the extinction estimates, especially if only a few 
stars are seen in the middle of the filament. We test the effect of 
YSOs by comparing results obtained before and after removing 
them. 

We also study t wo other filaments m apped with Herschel, 
already presented in lJuvela et al.1 (l2012bl) . The profiles of these 
filaments are shown in AppendixlAl 

2.1. WFCAM 
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Fig. 1. Extinction maps of TMC-1 (southern field, 2MASS data) 
and TMC-1N (northern field, WFCAM data). 



sample a 1° x 1° field in the near-infrared (NIR) J-, H-, and K- 
bands (1.2, 1.6, and 2.2 ftm, respectively) in the Taurus molec- 
ular cloud complex field TMC- IN. Our target field choice was 
based on the Tauru s extinction maps of lCambresvl (fl999) and 
IPadoan et all d2002h . 

We used the ON-OFF method in order to be able to measure 
diffuse emission. The observations were made during 13 nights 
between 2006-2008 on 4 x 4 adjacent ON-fields. 

The integration times of the ON-fields were 8.5, 11.3, and 
27.3 minutes in the J, H, and K band, respectively. The coordi- 
nates of the 4 OFF-fields are (RA (J2000), Dec (J2000)): 3h52m 
and +23°43', 3h52m and +24°9', 3h54m and +23°43', 3h54m 
and +24°9'. 

The data were reduced in the normal pipeline routineQ. In or- 
der to study the surface brightness, standard decurtaining meth- 
ods were needed to remove the stripes caused by the instrument. 
The reduced data wer e converted to 2MASS magnitude scale 
using the formulas of iHewett et al.l (|2006) and calibrated us- 
ing 2MASS stars. The Iraf Daophot-package was used for pho- 
tometry and Montag^l for combining the images. The median 
magnitudes for which magnitude uncertainty cr ~ 0.1"' of our 
WFCAM data are J = 20.0 m , H = 19. l m , and K = 18.8 m . These 
are 4 - 4.5™ deeper than 2MASS observations. In the analysis, 
we used only sources with cr < 0.1 m . Saturated stars were re- 
placed with corres ponding 2MASS sta rs. Based on the classi- 
fication scheme of ICantiello et alJ d2005h . potential background 
galaxi es were removed from our star list, in a manner similar to 
that of lSchmalzl eTaT] d2010h . 

2.2. 2MASS 

We used Two Micron All Sky Survey (2MASS) (Skrutskie et al. 
2006) data to make extinction maps of the fields for which we 
did not have deeper observations. We also used 2MASS data for 
TMC- IN in order to perform a comparison with the WFCAM 



We h ave used the Wide Field CAMera (WFCAM) dCasali et alj 
120071) of the United Kingdom InfraRed Telescope (UKIRT) to 



1 http : //casu. ast . cam. ac .uk/surveys-projects/wfcam/technical 

2 http : //montage . ipac . caltech . edu/ 
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data that contains a significantly larger number of stars. The sizes 
of our fields are approximately 1 x 1 ° . The limiting magnitudes 
of 2MASS are J = 15. 8 m , H = 15. l m , and K s = 14.3 m . 

2.3. Herschel 

Both TM C-1 and TMC-1N have also been mapped with 
Hersch el dPilbratt et all 1 201 Ol) as part o f the Gould Belt 
Survey (lAndre et al.ll201ol) . We used SPIRE (iGriffin et alJ bOlO) 
250//m, 350/mi, and 500 /vm maps of both of the fields. We ob- 
tained the data from the Herschel Gould Belt Survey consor- 
tium. The observation identifiers of the data are 1342202252 and 
1342202253. 

Fields G300.86-9.00 (PCC550) and G163.82-8.44have been 
mapped with Herschel as part of the Cold Core s project and are 
described in more detail in lJuvela et al.l (1201 2bl) . 

3. The column density maps 

3.1. Column densities from the NIR reddening 

The extinction maps were derived from the near-infrared red- 
dening of the background stars, using 2MASS observations 
for TMC-1 and the WFCAM observations for the northern 
field TMC-1N. For comparison, an extinction map of TMC- 
1N was also made using 2MASS data. Th e extinction values 
were calculated with the NICER method (Lombar di & Alvesl 
1200 assuming an extinction curve corresponding to R\ = 4.0 
(ICardelli et alll 19891) . 

For comparison with the N(H2) values below, the extinc- 
tion values were conv erted to N(H?), a ssuming a relation 
N(U 2 ) ~0.94xl0 21 A v dBohlin etalJll978h . The relation itself 
is derived for diffuse regions, but serves as a good point of refer- 
ence for the comparison between the measurements of dust red- 
dening and dust emission. The resolution of the derived column 
density map is 2' for 2MASS data and 1' for WFCAM data. 

3.2. Column densities from dust emission 

The maps of the dust colour temperature, 7c, were calculated 
using the SPIRE 250//m, 350//m, and 500 yum maps. The data 
were convolved to a common resolution of 40" and the Tq val- 
ues were estimated pixel by pixel assuming a constant value of 
the dust spectral index, B = 2.0. The least-squares fit of modified 
black body curves B v (7 , c)v e provides both the colour tempera- 
ture and the intensity that were converted to estimates of column 
density. The column density was obtained from the formula 



JV(H 2 ) 



B y (T)K/im H 



(1) 



where the dust opacity per total mass, k, was assumed to fol- 
low the law 0.1 cm 2 /g (v/lOO OGHz)^ that shou l d be applicable 
to hig h density environments (Hildebrand 1983; Bec kwith et alj 
1990). The resolution of the column density maps derived from 
Herschel observations is 40". 



4. Modelling the filament profiles 

We selected the most prominent filament from each field. The 
filaments were originally traced by eye and the paths of the fila- 
ments were consequently adjusted to follow the crest of the col- 
umn density as derived from the observations of the dust emis- 
sion. The filament profiles perpendicular to the filament are ex- 



tracted at small intervals and the final profile of the filament is 
obtained by averaging all data along the filament. 

We fit the filament column density profiles N(r) with 
Plummer-like profiles 



P P (r) = 



N(r) = A p 



PcRnat 



[1+Wfl flat ) 2 p- 1)/2 



(2) 



(see e.g. lNutter et al.l2008HArzoumanian et al .11201 lh . The equa- 
tion includes the central density p c , the size of the flat inner part 
Rf\ dt , and the parameter p that describes the steepness of the pro- 
file in the outer parts. We include in the fits two addit ional pa- 
ramete rs to describe a linear background. According to Ostriker 
the value of p should be 4 in the case of an isothermal 
cylinder in hydrostatic equilibrium. However, the values that are 
derived fro m observations are often smaller and typically a round 
p ~2 (e.g. [Arzoumania n et al. l l20Tll Uuvela et al.|[2012bb . The 



factor A p is obtained from the formula A p 



1 r°° 

COS! J— CO (1 



(III 



+U 2 )P' 2 ' 



where we assume an inclination angle of i - 0. In addition to the 
fitted Plummer profile parameters, we also derive the mass per 
unit length Mi; ne of the filament by integrating column density 
over the profile. 

We use four different methods for extracting the profiles, i.e., 
profiling. The methods are described in more detail in the follow- 
ing subsections. In short, the profiles are based on a fit to 

1 . a column density map derived from (submillimetre) dust 
emission map [Method A]. 

2. a column density map derived from NIR extinction map 
[Method B]. 

3. median Ay of stars within a certain offset from the filament 
centre [Method C]. 

4. A v of individual stars [Method D]. 

Method A is currently often used for estimating filament prop- 
erties. The aim of our study is to investigate, if the alternative 
methods B, C, and D based on NIR data could also be used for 
extracting filament profiles. 

4.1. Extraction of profiles from column density maps 

Methods A and B both use a column density map for fitting, but 
in method A, the column density map is derived from a (sub- 
millimetre) dust emission map, and in method B, from a NIR 
extinction map. Otherwise, the process for extracting the profile 
parameters is the same. We step along the filament with a step 
size equal to half of the full width at half maximum (FWHM) 
of the input map. At each position, we construct the local pro- 
file that consists of the values read from the map at intervals of 
0.3xFWHM along a line perpendicular to the filament length. 
Before averaging, the profiles are aligned in order that their 
peaks match on the offset axis. 

Because the observations have a finite resolution, the data are 
fitted with a Plummer profile that is convolved to the resolution 
of the column density map used. Thus, the obtained parameters 
of the Plummer function should correspond to those of the true 
filament. In addition to the convolved Plummer profile, we in- 
clude in the fits a linear background to take into account the pos- 
sible extended background. The final free parameter is the shift 
of the Plummer profile that allows for some errors in the posi- 
tion of the originally selected filament. For the dust emission, 
this already follows the maximum column density. However, it 
is better to allow a shift in the fit. This applies particularly to the 
analysis of the extinction where, in principle, the exact position 
of the observed filament could differ from the position deduced 
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from the dust emission (as is the case in TMC-1, see Figs.|6]and 
0. 

The estimated FWHM of the filament profile and mass per 
unit length, My me , depend on how wide a region is used in the 
fitting and how the background is removed, as there are no clear 
limits of where the filaments end, particularly in crowded ar- 
eas. The fitted Plummer parameters, however, should not depend 
on the width of the fitting region, as long as the region is wide 
enough to cover also the edges of the filament. In the fits to the 
average profile, we only use data where offsets from the fila- 
ment centre correspond to distances [-0.3 pc, +0.3 pc] for TMC- 
1 and TMC-1N. For the other filaments shown in Appendix lAl 
we use bigger offset limits, 0.5 pc, in order to make sure that 
these wider and more distant filaments are properly fitted. We 
calculate FWHM values for each individual profile and for the 
median profile after subtracting the minimum. This is directly 
the observed value without deconvolution with the instrument 
beam. 



4.2. Extraction of profiles using Ay of stars directly 

The previous extinction maps represent the extinction values 
of the individual stars that are averaged with a relatively large 
Gaussian beam. This is necessary in order that no holes remain 
in the Ay map and each pixel value corresponds to the average 
of at least of a few stars. However, because we are mainly inter- 
ested in the average properties of the filaments, we can average 
all the data along the filament length and this should allow us 
to reach a higher resolution for the profile function. We use two 
different methods for extracting profiles: median Ay values for 
each offset bin (method C) or direct fit to the Ay values of in- 
dividual stars (method D). In both of these methods we extract 
profiles only using stars that are within a given offset from the 
filament centre: [-0.3 pc, +0.3 pc] for TMC-1 and TMC-1N, and 
[-0.5 pc, +0.5 pc] for the other filaments shown in Appendix lAl 
as in the case of using column density maps. 

In method C, we calculate the offset from the filament cen- 
tre for each star. The stars are binned in narrow offset intervals 
(i.e., bins extended along the whole filament) and we calculate 
the median Ay for each offset with a step size 0.3 times the bin 
width. These profiles can again be fitted, taking into account that 
the binning effectively introduces some loss of resolution. In the 
following results and simulations, the width of the offset bins 
is always one arc minute, although we have tested also other 
values. The drawback of this method is that we cannot use the 
average and the dispersion of the Ay values of the neighbouring 
stars to filter outliers as effectively. We calculate the FWHM of 
the median profile after subtracting the minimum. As with col- 
umn density maps, this is directly the observed value without 
deconvolution with the instrument beam. 

Method D relies directly on the measurements of the indi- 
vidual stars, without any spatial averaging. Each star represents 
extinction within a very narrow region and the final resolution of 
the profile function depends only on the density of the stars along 
the offset axis. The influence of each star is weighted not only 
according to the uncertainty of its photometry but also inversely 
proportionally to the local stellar density around the star. Thus, 
the final profile is weighted according to the area rather than the 
number of stars within that area. The data are fitted directly with 
the Plummer profile. Because the individual stars provide a prac- 
tically infinite resolution, no convolution is needed. In connec- 
tion with this method we do not calculate any separate FWHM 
estimates. 



5. Results 

5. 1 . Analysis of the column density maps 

Column density maps and filament profiles derived from the 
WFCAM extinction map and Herschel emission map of TMC- 
1N are shown in Figs. [2]and[3] The filament is traced based on 
the Herschel map and marked with white line in both column 
density maps. The filament position obtained from the Herschel 
map fits very well also in the WFCAM map. In the Herschel map 
the marked part forms the densest continuous filament, which we 
call here the full filament. However, in the WFCAM extinction 
map, the filament continues seamlessly southwards outside the 
frame. 

The profiles obtained by using these two methods (A and 
B) are very similar. The FWHM of the filament is approx- 
imately 0J__2C J _j^rMi_J^in_a£reement with the general re- 
sults of lArzoumanian et al.l d201 ll) on filaments in the IC 5146 
molecular cloud. When plotting together the column density and 
FWHM values along the ridge of the filament, we can see a 
notable anticorrelation between these parameters. In the dens- 
est parts of the filament FWHM gets smaller, possibly indi- 
cating that gravity is pulling together the mass and forming 
na rrower condensation s in the filament, similarly to the results 
of lJuvela et aT] d2012bl) . 

The column density structure along the ridge of the filament 
is quite similar in both data sets. In the Herschel data, the highest 
peak is narrower and reaches a column density of ~ 30 • 10 21 
cm' 2 while in the WFCAM data the maximum is ~ 25 ■ 10 21 
cm -2 . Considering the many uncertainties in the column density 
estimates (e.g., dust opacity), the values are practically identical. 
Conversely, the similarity of the column density values suggests 
that the assumptions used in the column density calculation were 
reasonable. 

To study the effect of choosing the filament end points, we 
also extract profiles from the middle part of the filament in TMC- 
1N, shown in Fig. |4] The average column density is higher than 
in the full filament. The fitted Plummer parameter values are 
slightly altered, although the results are quite close to the val- 
ues obtained using the full filament. 

For a comparison using lower resolution data, a column den- 
sity map and filament profile derived from the 2MASS extinction 
map of TMC-1N are shown in Fig. [5] Here, the filament position 
obtained from the Herschel map also fits the filament seen in 
2MASS extinction map very well. However, along the filament, 
the column density peaks are in different positions in maps de- 
rived from the Herschel and 2MASS observations. Possible anti- 
correlation between FWHM and column density along the ridge 
of the filament is not evident using the 2MASS data. The col- 
umn density derived from the 2MASS data is notably smaller 
than the column density derived from the Herschel or WFCAM 
data. However, we also obtain a similar Plummer profile using 
2MASS data. 

For a comparison of another filament, column density maps 
and filament profiles derived from the 2MASS extinction map 
(YSOs removed) and Herschel emission map of TMC-1 are 
shown in Figs. [6] and [7] Here the filament is again traced us- 
ing the Herschel column density map. In this case, the filament 
position does not exactly match the filament ridge in the 2MASS 
map. However, this is not a problem in the profiling process, as 
the filament is allowed to have a small offset from the given po- 
sition. In this filament, we do not see any notable background 
gradient or anticorrelation between FWHM and column density 
along the ridge of the filament. The FWHM of the filament is 
approximately 0.1-0.2 pc. The column density structure along 
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Fig. 2. TMC-1N column density map and filament profile derived from WFCAM extinction map. a) Column density map. The part 
of filament used in the fitting is marked with white line (based on the Herschel column density map in Fig. [3}. b) Average column 
density profile of the filament (black line), fitted Plummer profile (red dashed line) and the baseline of the fit (blue dashed line). 
Values for fitted Plummer parameters p c , Z?fl at , and p are marked in the figure, c) FWHM values (black circles and left hand scale) 
and column density along the ridge of the filament (solid line and the right hand scale). 
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Fig. 3. TMC-1N column density map and filament profile derived from Herschel emission map. See Fig.|2]for explanations of the 
used notation. 



the ridge differs notably in these two data. In the Herschel data, 
we see a continuous high density ridge, while in the 2MASS 
data there is more variation, because of the low number of back- 
ground stars. Based on the Herschel data, next to the filament 
there is a distinct clump, which effects the profile and the de- 
rived FWHM values. This change in the FWHM values is not 
seen in the 2MASS data. For a comparison, we show the col- 
umn density map and filament profiles derived from the 2MASS 
extinction map before removing YSOs in Fig. IB. II 



The obtained parameters describing the filament, namely 
FWHM, Plummer parameters and the mass per unit length Mii ne 
derived from observations of TMC- 1 N (full and middle filament) 
and TMC-1 (before and after removing YSOs) using different 
methods are shown in Table Q] 

We can compare our results for TMC-1 deriv e d from 
the Herschel maps with those of Nutte r et al ] d2008), 
from Submil limetre Common User Bolometer Array 
(SCUBA) dHolland et alJ [[999) and the Infrared Space 
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Fig. 4. TMC-1N column density map and profile of the densest middle part of the filament derived from WFCAM extinction map. 
See Fig.[2]for explanations of the used notation. 
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Fig. 5. TMC-1N column density map and filament profile derived from 2MASS extinction map. See Fig.|2]for explanations of the 
used notation. 



Obser vatory ISO dKessler et all 1 19961) ISOPHOT instru- 
ment dLemke et al.lll996l) data. We note first that we agree on 
the fact that the Bull's tail, as they call it, is a filament and not a 
ring. Both analyses fit Plummer-like profiles, and therefore the 
results should be comparable (see Table [TJ. We obtain a value 
for p of ~2.3, whereas they fix p=3. We find p c ~ 7 x 10 4 , 
compared to their range of p c ~ 7 x 10 4 to ~ 1.7 x 10 5 (with 
nominal values towards the upper end of this range). Our value 
for Rfia, is ~ 0.04pc, and their value is ~ 0.05pc. We find the 
FWHM to be 0.1 2pc for the Herschel data, and -0.1 8pc for the 
2MASS data. They fitted the filament with two components, the 
first of which had a FWHM of 0.1 lpc and the second of which 



had a FWHM of 0.16pc. It can therefore be seen that there is a 
high degree of consistency, giving us added confidence in our 
parameter fitting of the data. 

The WFCAM Ay maps and Herschel column density maps 
give very similar results for the FWHM of the TMC- 1 N filament 
(0.092 and 0.093 pc, respectively). Lower resolution Ay maps 
based on 2MASS stars give slightly higher values for FWHM 
(0. 17 pc). The behaviour of the Plummer parameters or My me de- 
rived from these is not as predictable, as there are dependencies 
between the parameters. We examine these dependencies more 
closely below. In addition, the exact values seem to be quite sen- 
sitive for instance to small changes in the data, different reso- 
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Fig. 6. TMC-1 column density map and filament profile derived from 2MASS extinction map (YSOs removed). See Fig. [2] for 
explanations of the used notation. 
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Fig. 7. TMC-1 column density map and filament profile derived from Herschel dust emission map. See Fig. [2] for explanations of 
the used notation. 



lution data, fitting algorithm, initial fitting parameter values or 
different profiling method. We also noticed that the chosen pixel 
size for the extinction maps can alter the results. 

In order to study the correlations of the Plummer profile pa- 
rameters we save cross-sections of x 2 - We fix each of the param- 
eters p c , /?fl at or p in turn to the fitted value, and change the others 
systematically to extract the;^ 2 values in three orthogonal param- 
eter planes crossing at the location of the x 2 minimum, Xmm- ^ n 
Fig.[8l we show the contour levels 1.5 y 2 and 10 y 2 , for each 
case. The results show that there is a banana shaped correlation 
between and p, more or less linear (in log-space) correlation 
between p c and p, and a linear anticorrelation between p c and 



Rf\. dt . Even the 1-5^- contours are quite wide, especially in p c - 
/?fl at -plane, meaning that very different parameter values can give 
small x 1 values. There are also differences in the size and orien- 
tation of the contours when comparing Herschel emission maps 
and WFCAM extinction maps. The l-5^ 2 ;n contours are smaller 
when using WFCAM extinction map in /?fl at -p-plane and p c -p- 
plane. However, in p c -/?fl at -plane the contours are smaller when 
using Herschel emission maps. 
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Table 1. Obtained parameters describing the filaments TMC-1N (full and middle filament) and TMC-1. The columns are the data 
set, method (column density map (A), Ay map (B), median profile (C) and fit to individual stars (D)), FWHM (including the effect 
of the data resolution), Plummer parameters and mass per unit length M]; ne . 



Data set 


Method 


FWHM 


Pc 




P 


Mine 






(pc) 


(cm- 3 ) 


(pc) 




(M Q /pc) 


TMC-1N: full filament 


Herschel maps 


A 


0.093 


3.94e+04 


0.035 


2.81 


1.91e+01 


WFCAM: A v maps 


B 


0.092 


9.02e+04 


0.012 


1.84 


3.04e+01 


2MASS: Am maps 


B 


0.17 


4.28e+04 


0.016 


1.67 


3.40e+01 


WFCAM: Median profile 


C 


0.12 


1.70e+04 


0.074 


4.57 


1.30e+01 


WFCAM: Fit to individual stars 


D 




2.17e+04 


0.058 


3.71 


1.49e+01 


2MASS: Median profile 


C 


0.16 


1.55e+04 


0.048 


2.44 


1.98e+01 


2MASS: Fit to individual stars 


D 




1.62e+04 


0.068 


3.76 


1.49e+01 


TMC-1N: middle filament 


Herschel maps 


A 


0.088 


5.49e+04 


0.028 


2.37 


2.93e+01 


WFCAM: A v maps 


B 


0.088 


1.67e+05 


0.007 


1.56 


5.95e+01 


2MASS: A v maps 


B 


0.21 


1.78e+04 


0.044 


1.93 


3.97e+01 


WFCAM: Median profile 


C 


0.15 


1.70e+04 


0.063 


2.92 


2.33e+01 


WFCAM: Fit to individual stars 


D 




4.33e+04 


0.023 


1.79 


4.40e+01 


2MASS: Median profile 


C 


0.18 


1.35e+04 


0.087 


3.38 


2.49e+01 


2MASS: Fit to individual stars 


D 




2.66e+04 


0.028 


1.63 


5.55e+01 


TMC-1 


Herschel maps 


A 


0.12 


6.98e+04 


0.043 


2.27 


9.02e+01 


2MASS: A v maps (with YSOs) 


B 


0.17 


3.91e+04 


0.078 


3.60 


5.06e+01 


2MASS: A v maps (without YSOs) 


B 


0.18 


3.73e+04 


0.062 


2.85 


5.24e+01 


2MASS: Median profile 


C 


0.05 


9.60e+04 


0.019 


1.84 


6.56e+01 


2MASS: Fit to individual stars 


D 




6.56e+04 


0.032 


1.95 


8.23e+01 



5.2. Profiles using Ay of stars directly 

Profiles of TMC-1N full filament derived by using methods C 
and D with WFCAM data and 2MASS data are shown in Fig. [9] 
The profile of TMC-1 N middle filament using the WFCAM data 
is shown in Fig.[l0] and the profile of TMC-1 using the 2MASS 
data is shown in Fig. Q~T] The obtained parameters derived from 
observations of TMC-1N (full and middle filament, WFCAM 
and 2MASS data) and TMC-1 (2MASS data) are shown in 
Table Q] The results of using median averages (method C) are 
on the 'Median profile' and the results of using individual stars 
(method D) on the 'Fit to individual stars' line. The results ob- 
tained for TMC-1 before and after removing the YSOs were 
identical when using methods C and D. 

The values obtained with methods C and D are quite sim- 
ilar to the values derived from maps. However, as with the 
map method, the derived parameters are quite sensitive to small 
changes in the data. For example, for the TMC-1N full filament 
using the WFCAM data, the star methods give notably larger val- 



ues for p than the map methods, but when only the densest mid- 
dle part of the filament is profiled, these differences get smaller. 

6. Simulated extinction observations 

We use simulations to examine the accuracy with which the fil- 
ament profiles can be estimated from the observations of back- 
ground stars. As a reference, we use the WFCAM observations 
from which we estimate the probability distribution of the num- 
ber of stars as a function of the H-band magnitude, the average 
J-H and H-K colours, and the uncertainties of these colours as 
a function of the magnitude. These data are used in simulations 
in which we place background stars at random positions over an 
area of 30'x30'. As a reference, we use the stellar density in the 
WFCAM observations but also test the effect of lower and higher 
stellar densities. The stellar colours are reddened according to 
the column density distribution that results from the presence of 
a filament. The density distribution of the filament follows the 
Plummer function with chosen parameters. The filament is as- 
sumed to be in the plane of the sky but has a position angle of 



8 



J. Malinen et al.: Profiling filaments: comparing near- 

Offset (pc) 



.-- p-tpf/pp-pp.- 




6 a 



--0 2 

Offset (arcmin) 
Offset (pc) 



0.0 



0.1 



0.2 




8 

-7-p- 
6 

i — - 
5 o 



3 ffi 
2 



- 2 

Offset (arcmin) 
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frame). Blue stars: median values. Gray line: Plummer profile 
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solid line: The convolved Plummer profile. Red solid line: fit to 
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30° in order to avoid the special case where the filament would 
be aligned with the axes used in the pixelization of the extinction 
maps. 

The extinctions of the individual stars (and of the 30'x30' 
extinction map) are calculated using one corner (78 square arc 
minutes) of the simulated region as a reference area from which 
the intrinsic stellar colours are estimated. The simulated obser- 
vations are then analysed as in the case of the WFCAM obser- 
vations above. The filament parameters are derived from the ex- 
tinction map (method B), from the median profile along the fila- 
ment (method C), and from the direct fit of the Plummer profile 
to the Ay values of the individual stars (method D). We anal- 
yse a 20' long part of the filament that thus fits well within the 
30'x30' area with the assumed distance of 140 pc, even taking 
into account the non-zero position angle. 

In the simulations, R^ t is always 0.025, and there is no back- 
ground gradient. We compare our three methods and the way the 



fitted p, p c and derived Mii ne change when we use different sets 
of parameter values for p c (10 4 , 10 5 , 10 6 cm -3 ), p (2, 3, 4), and 
stellar density (0.2, 0.5, 1, 2, 5, 10, where 1 denotes the stellar 
density in the WFCAM observations). We also test how photo- 
metric noise affects the results by using different noise levels 0, 
0.3, and 1, where 1 denotes the photometric noise level in the 
WFCAM observations of TMC-1N. Even when the photometric 
noise is zero, there is still noise in the data, because the stars 
have a random position and brightness, i.e., some stars may be 
extincted below the detection threshold. The filament looks dif- 
ferent depending on the number and brightness of the stars sit- 
uated in the filament area. The filament is well defined, if there 
are enough bright stars near the borders of the filament, where 
they can still be seen. We run each simulation case 100 times for 
the estimation of the uncertainties. 
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6.1. Simulation results 

Examples of simulated maps and filament profiles with and with- 
out photometric noise are shown in Figs. [T2land[T3l The sim- 
ulated profiles of the same cases using methods C and D are 
shown in Fig. [14] The resulting column density maps resem- 
ble our TMC-1N observations, except that the density distribu- 
tion precisely follows the Plummer function with a chosen set 
of parameters. If all the simulated filaments are examined along 
their whole length, there is no clear anti-correlation (nor cor- 
relation) between observed FWHM and column density along 
the ridge. However, similarly to the observations of TMC-1N 
(Fig. [2] frame c), anticorrelation can sometimes be seen in some 
of the densest peaks also in simulations, even though the simu- 
lated filament itself is symmetric. 

The filament parameters are derived from the simulated ob- 
servations. In Figs. 15-20, we compare the true and estimated 
parameter values p, p c , and Mr me for cases with different param- 
eter values and different number of background stars. The error 
bars correspond to half of the interquartile range between the 
75% and 25% percentiles. 

Figure[l5]shows that in a noiseless situation method D works 
almost perfectly in fitting the p parameter. The error bars of 
method D are not visible in the figure, because they are so small. 
This is natural as the random positions of the stars do not af- 
fect this method and the photometry itself does not contain er- 
rors. Methods B and C also work well for modest density, p c = 
10 4 cirT 3 , being correct to within 10%. However, in higher den- 
sity (p c = 10 5 cm 4 ) both of these methods underestimate the p 
value, except in the case p — 2, where method B ends up over- 
estimating the p value. Still, both methods are correct to within 
20%. Method B does not work reliably at densities as high as 
10 6 cm 4 , because stars cannot be seen at the centre of the fila- 
ment. It is therefore left out of the figures. The case of density 
10 6 cm 4 and p — 2 is not shown either with methods C and 
D, as the results get unreliable when the empty part in the mid- 
dle is too wide. However, with larger p values both methods C 
and D continue to work also at high densities, although method 
C underestimates the p value ~20%. We also made some tests 
with density 10 3 cm 4 , but for such low density the filament is 
no longer reliably identified. The combination of central density 
and true parameter p determines how well the observed param- 
eter p is fitted. 

The results with photometric noise similar to WFCAM ob- 
servations are shown in Fig. Q~6] The results are still similar to 
those without noise, except that the errors get larger. For density 
10 4 cm 4 and parameter p-2.0 the bias is very small with all star 
densities and all methods, but in the case of very low star density 
(0.2), the scatter for method C is so large that the method does 
not work reliably. For density 10 4 cm 4 and parameter p=4.0, 
both methods B and C typically result in a value of p = 10, the 
maximum allowed in the fitting procedure. This means that in 
these cases the results of the fits are not reliable. Still, in most 
cases the results are correct to within 10% for method D and to 
within 20% for the other methods. 

Figures [17] and [18] show, that p c is also fitted quite reliably, 
especially using methods D and C for modest (p c ~ 10 4 cm 4 ) 
densities and without photometric noise. If noise is included, 
method B tends to overestimate the density for the modest den- 
sity case, especially when the star density is low. For higher den- 
sities (p c >= 10 5 cm 4 ), both methods B and C underestimate the 
density, even though the results are still correct to within an or- 
der of magnitude. The derived Mi; ne shows a similar behaviour in 
Figs.[19]and [20] Method D gives correct results to within 10% in 



most cases. Methods B and C also give correct results to within 
10% for modest density (p c ~ 10 4 cm 4 ). For higher densities, 
methods B and C both underestimate the mass by ~20% or even 
up to ~60% in the densest case. 

Generally, all the compared methods seem to produce reli- 
able results in most cases, especially with modest density (10 4 
cm 4 ). In addition, when stellar density grows, the errors be- 
come smaller. For higher densities, method D usually seems to 
give the most reliable value. 
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Fig. 12. Example of a simulated map and filament profile without noise. Used parameter values in the simulation: p c =10 4 cm~ 3 , 
^flat=0.025, p-2.Q, stellar density = 1.0, and without photometric noise. The fitted parameter values are shown in frame b. See 
Fig.[2]for explanations of the used notation. 
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Fig. 13. Example of a simulated map and filament profile with noise. Used parameter values: p c = 10 4 cirT 3 , /?fl at =0.025, p-2.0, stellar 
density = 1.0, and noise=1.0 (meaning a similar noise level as in our WFCAM observations of TMC-1N.). The fitted parameter 
values are shown in frame b. See Fig. |2]for explanations of the used notation. 
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Fig. 14. Examples of simulated profiles using methods C and D. Used parameter values: p c =10 4 ctrr 3 , /?fl at =0.025, p=2.0, stellar 
density = 1.0, and without noise (left frame) or with noise (right frame). The fitted parameter values are shown in the figure. See 
Fig.[9]for explanations of the used notation. 
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7. Discussion 



The characteristics of interstellar filament properties are cur- 
rently approached using mainly submillimetre studies, due to the 
high resolution data obtained with Herschel. However, if such 
data are not available, the use of groundbased NIR data would 
be more convenient. We have therefore compared the profiles 
obtained from column density maps derived from both submil- 
limetre dust emission maps (method A) and NIR extinction maps 
(method B). 



We have presented two new methods to derive filament pro- 
files. In method C we divide the stars to bins according to their 
distance from the filament centre, calculate the median Ay for 
each offset bin and fit a Plummer profile to the median Ay val- 
ues. In method D we fit the Plummer profiles directly to the Ay 
values of individual stars, without any spatial averaging. 

The results of the profile fits to TMC-1N and TMC-1 
(Table [TJ show, that all four methods give rather similar re- 
sults. Of course, with observations we do not know which val- 
ues are the most correct ones. There does not seem to be any 
obvious systematic biases in different methods, e.g., Ay maps 
can sometimes give a lower (in TMC-1) and sometimes higher 
(in TMC-1N) values for masses than the Herschel emission 
maps. In general, the fitted densities are similar to within an or- 
der of magnitude. The fitted p values change relatively much, 
between 1.67-4.57 for TMC-1N (full filament), 1.56-3.38 for 
TMC-1N (middle filament), and 1.84-3.60 for TMC-1, depend- 
ing of the method. However, in the other two filaments shown in 
Appendix lAl the p values are much more stable (1.92-2.1 and 
1.31-1.79). 

YSOs embedded in the filaments could bias the extinction 
estimates, as they are redder than background stars. Even the 
effect of a single YSO could be notable in dense filaments, if 
only a few stars are seen in the middle of the filament. There 
are six known YSOs in the TMC-1 filament area used in our 
profiles. We have tested the effect of YSOs in this filament, and 
note that removing the YSOs can have significant local effects 
in the column density maps and profiles along the ridge of the 
filament. However, in this case, the profile parameters derived 
with or without YSOs stay identical when using methods C or 
D. Using method B, there is some change in the fits, but the 
parameter values are still similar to within ~ 10-30%. 

We have noticed that the Plummer profile fits are quite sen- 
sitive to small changes in the data or the fitting routines. The 
dependencies of the Plummer parameters in Fig.|8]show that the 
confidence regions of Plummer parameters are quite wide. This 
confirms our findings, that the parameters are not always well 
constrained and therefore small changes in the data can lead to 
different parameter values. The correlations of the parameters, 
e.g. 7?fl at and p, should therefore be taken into account, instead 
of giving too much weight on the value of a single parameter. 
The emission and extinction data are independent and together 
should result in more accurate estimates of the filament parame- 
ters. This is particularly true when the orientations of the corre- 
sponding confidence regions (i.e., the minimum^ 2 contours) are 
different. 

We have used simulations to compare the filament param- 
eters and the derived mass per unit length Mii ne obtained with 
extinction maps and the two alternative methods. We have tested 
the limit of the methods using different densities p c (10 4 , 10 s , 
10 6 cm 3 ) and stellar densities (0.2, 0.5, 1, 2, 5, 10, where 1 
denotes the stellar density in the WFCAM observations). Tests 
with p c = 10 3 showed, that in these low densities the extinc- 
tion methods did not work reliably. The simulations show that 
method D gives correct results to within 10% in most cases. 
Methods B and C also work in most cases correctly to within 
~20%. Only for high densities (10 6 ctrT 3 ) can the bias for 
method C be notably larger, and method B cannot reliably be 
used. When using method D in high density the scatter can grow 
so much that the results are unreliable, even though the bias 
can still be small. When stellar density grows, the results be- 
come more accurate. However, the simulations showed that all 
the methods also work in most cases with stellar densities below 
that of our WFCAM observations. 
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Fig. 15. Plummer p parameter derived from simulated extinction observations without noise. Different central densities are shown 
on different lines: 10 4 (cm -3 ) (uppermost frames), 10 5 (middle frames), and 10 6 (bottom frames). Symbols and colours denote the 
used method: method B (black square), method C {blue circle), and method D {red star). Stellar density is given on the x-axis, 
where 1 means a similar density as in the WFCAM observations of TMC-1N. The true p value (2, 3 or 4) used in the simulations is 
given on top of each column and also shown as a dashed line in the figures. The derived p value is marked in the figures for each 
simulated stellar density and method with scale on the y-axis. 
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Fig. 16. Plummer p parameter derived from simulated extinction observations with noise = 1 (times WFCAM noise). See Fig.[T31 
for explanations of the used notation. 
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Fig. 17. Plummer p c parameter derived from simulated extinction observations without noise. Y-axis denoting the fitted p c is divided 
with the correct p c . See Fig.[l5]for explanations of the used notation. 
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Fig. 18. Plummer p c parameter derived from simulated extinction observations with noise 
of the used notation. 
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Fig. 19. Mii n e derived from simulated extinction observations with noise = 0. The true Mi; ne is marked with a dashed line in the 
figures. The derived Mn ne is marked in the figures for each simulated stellar density and method with scale on the y-axis. Otherwise, 
see Fig.[l5]for explanations of the used notation. 
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Fig. 20. Mii n e derived from simulated extinction observations with noise = 1. See Figs. | 
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In nearby, high galactic latitude clouds, such as Taurus, fore- 
ground stars do not cause a problem in the analysis. In other 
cases, however, foreground stars could cause a notable bias. 
We have examined the effect of foreground stars in simulations. 
Changing ~10% of the background stars to foreground stars can 
cause ~ 10-30% more bias in the mass estimates. The method D 
is quite robust when no noise is included, but the difference to 
the other methods is diminished, when realistic noise is added. 

The key question is, what kind of conclusions can be made 
of the observations on the basis of simulations. In the WFCAM 
observations of TMC-1N, the value of central density is esti- 
mated to be between 10 4 -10 5 cirT 3 , which is within the range of 
our simulations. In this case, in particular method D should be 
quite reliable according to our simulations (see Figs.[T6land l20b . 
However, the simulations are simplifications of true filaments, 
and they do not take into account all factors that could affect 
the interpretations. For example, there is no background gradi- 
ent in the simulations, whereas in the observed area, there is a 
strong gradient that could distort the obtained parameters. A lin- 
ear background is included as a free parameter in the fitting of 
Plummer profile, and therefore a linear background gradient per- 
pendicular to the filament should not cause problems. However, 
nonlinear gradients could affect the fitting. In the simulations the 
filament is also assumed to be in the plane of the sky. All the gen- 
erated photometric errors are normally distributed. In real obser- 
vations, there could be outliers, which could affect both methods 
C and D. 

However, the fundamental question is, whether or not the 
real filament profiles are really Plummer functions. In the simu- 
lations, we presume that the filament profile obeys the Plummer 
function, and therefore the results fit well that function unless 
there is some bias. In real observations, this is not necessarily 
the case. We have also noticed that even a poor fit can give rea- 
sonable values for the Plummer parameters. Conversely, good 
fits can also give, e.g., very high values for the p parameter. It is 
therefore advisable that all fits are checked by eye, even for large 
surveys. 

The methods could be further improved, e.g., by carrying out 
a separate outlier filtering step (on Ay map) before using method 
D, or weighting the profiles according to the number of stars 
in the median bins in method C. Background structure could 
also be modelled and subtracted before extracting the profiles. 
In principle, a linear background should not affect the resulting 
profiles, but in practise the fits seemed to be quite sensitive to 
any small variations in the data. 

We have shown that filament profiles can be derived from 
NIR extinction data. These methods can serve as a valuable point 
of comparison to the results obtained with Herschel or other in- 
struments using submillimetre dust emission. The methods can 
be used with, e.g., UKIDSS and future NIR surveys, but we have 
also shown that even 2MASS data can give rather reliable es- 
timates in the case of very nearby clouds. NIR extinction and 
groundbased dust emission measurements with, e.g., SCUBA-II 
can be used as alternative or complementary means of study- 
ing the properties of filaments. The limits of the NIR extinction 
methods come from the requirements of having enough back- 
ground stars. This restricts the use of these methods to low lati- 
tudes and nearby sources. The resolution of the method is limited 
only by the amount of background stars. 

In order to achieve b etter resolution, NIR surface brightness 
could be used (see e.g. Pad oan et al.l l2006t lJuvela et all 120061 
l2008h . However, the observations of NIR surface brightness can 
be quite difficult. We study the properties of TMC-1 N using NIR 



surface brightness data in a following paper (Malinen et al. 2012, 
in prep.). 

8. Conclusions 

We have compared the filament profiles obtained by submillime- 
tre observations (method A) and NIR extinction (methods B, C, 
and D). 

- We conclude that NIR extinction maps (method B) can be 
used as an alternative method to submillimetre observations 
to profile molecular cloud filaments. 

- We show that if better resolution NIR data are not available, 
2MASS data can also be used to estimate approximate pro- 
files of nearby filaments. 

- We also present two new methods for the derivation of fila- 
ment profiles using NIR extinction data: 

• Plummer profile fits to median Ay values of stars within 
certain offset bin from the filament centre (method C) 

• Plummer profile fits directly to the Ay values of individ- 
ual stars (method D) 

- We show that in simulations, all the methods based on NIR 
extinction work best with modest densities p c = 10 4 — 10 5 
cm -3 . Method D gives correct results for the fitted and de- 
rived parameters to within ~10% and methods B and C to 
within ~20% in most cases. At high densities (~ 10 6 cm -3 ), 
only method D continues to work reliably. 

- For the profile fits to real observations the values of individ- 
ual Plummer parameters are in general similar to within a 
factor of ~2 (in some cases up to a factor of ~5). Although 
the Plummer parameter values can show significant varia- 
tion, the derived estimates of filament mass (per unit length) 
usually remain accurate to within some tens of per c ent. 

- We confirm the earlier results of N utter et al.l (t2008b that the 
examined TMC-1 structure, or Bull's tail, is a filament and 
not part of a continuous ring. Our parameter values are also 
consistent with the earlier results. 

- The confidence regions of the Plummer parameters are quite 
wide. As the parameters are not always well constrained, 
even small changes in the data can lead to different parameter 
values. The correlations of the parameters should therefore 
be taken into account. 

- The size and orientation of the confidence regions are not 
identical when Herschel maps or NIR data are used. The 
combination of the methods can significantly improve the 
accuracy of the filament parameter estimates. 
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Fig. A.5. Profiles of 300.86-9.00 (PCC550) filament derived by 
using methods C and D on 2MASS data. See Fig.[9]for explana- 
tions of the used notation. 
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Fig.A.6. Profiles of G163.82-8.44 filament derived by using 
methods C and D on 2MASS data. See Fig. [9] for explanations 
of the used notation. 



Appendix A: Comparing filament profiling in 
Herschel fields and 2MASS data 

We compare our methods to two profiles o btained us- 
ing He rschel column density maps presented in lJuvela et"aTI 
d2012bl) . namely fields G300.86-9.00 (PCC550) and G163.82- 
8.44. The distances to these filaments are ~230pc and 350 pc, 
respectively. Column density maps and filament profiles derived 
from 2MASS extinction maps and Herschel emission maps are 
shown in Figs.[AUand [O|for G300.86-9.00, and|A3]and|Al 
for G 163. 82-8. 44. 

The filaments are traced based on the Herschel maps and 
marked with white line in both column density maps. The fila- 
ment positions obtained from Herschel maps fit quite well also in 
the 2MASS maps. Profiles of the filaments derived by methods 
C and D using 2MASS data are shown in Figs. IA.5l and lA.6l The 
filament parameters derived with different methods are shown in 
table I A. ll The differences to the results of lJuvela et al.l 
are caused by small modifications in the method and the change 
of the selected filament segment. The distance to these filaments 
is greater than to TMC-1, which means that there are fewer stars 
in the filament area and a larger probability of foreground stars, 
which can bias the results. However, in general, we obtain quite 
similar results using all our four methods. In these filaments an- 
ticorrelation of FWHM in the densest peaks is also quite clear in 
the Herschel maps, but not in the 2MASS maps. 



Appendix B: Effects of YSOs 

YSOs may bias the extinction estimates, if they are not removed. 
For a comparison to Fig. |6]made after removing the YSOs (six 
in the profiled filament area), we show in Fig. IB. II the column 
density map and profiles of TMC-1 made before removing the 
YSOs. Even though the derived profile parameter values for both 
cases are similar to within ~10-30%, notable local changes can 
occur in the column density maps and in the column density pro- 
files along the ridge of the filament. 
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Fig.A.l. G300.86-9.00 (PCC550) column density map and filament profile derived from 2MASS extinction map. See Fig.|2]for 
explanations of the used notation. 
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Fig.A.2. 300.86-9.00 (PCC550) column density map and filament profile derived from Herschel emission map. See Fig. [2] for 
explanations of the used notation. 
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Fig. A.3. G163. 82-8.44 column density map and filament profile derived from 2MASS extinction map. See Fig.|2]for explanations 
of the used notation. 




N(H 2 ) (10 21 cm— 

1.5 3.0 4.5 6.0 7.5 9.0 10.5 12.0 
■ i i r 



in 



Offset (pc) 
-p-ip-p 77p-°p0.00 0.15 0.30 0.45 



24'00.0" 



12'00.0" 



o +37 °pppp- 
o 



48'00.0" 



I 

a. 



+36° pp-p 




-t -- 2 4 

Offset (arcmin) 
5 10 15 20 25 30 35 40 




30m24.00s 28m48.00s 4h27ml2.00s 
RA (J2000) 



0.0 1= i i i i i i i i I 2 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 

Offset (pc) 



Fig. A.4. G163. 82-8.44 column density map and filament profile derived from Herschel emission map. See Fig.|2]for explanations 
of the used notation. 
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Table A.l. Obtained parameters describing the filaments G300.86-9.00 and G163.82-8.44. The columns are the data set, method 
(column density map (A), Ay map (B), median profile (C) and fit to individual stars (£>)), FWHM (including the effect of the data 
resolution), Plummer parameters and mass per unit length M\ me . 
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Fig. B.l. TMC-1 column density map and filament profile derived from 2MASS extinction map (YSOs not removed). See Fig.|2]for 
explanations of the used notation. See Fig. |6]for a map and profiles derived after the YSOs are removed. 



